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ABSTRACT 

We study density cusps in the center of clusters of galaxies to reconcile X-ray mass 
estimates with gravitational lensing masses. For various mass density models with cusps we 
compute X-ray surface brightness distribution, and £t them to observations to measure the 
range of parameters in the density models. The Einstein radii estimated from these density 
models are compared with Einstein radii derived from the observed arcs for Abell 2163, Abell 
2218, and RX J1347.5-1145. The X-ray masses and lensing masses corresponding to these 
Einstein radii are also compared. While steeper cusps give smaller ratios of lensing mass to 
X-ray mass, the X-ray surface brightnesses estimated from flatter cusps are better fits to the 
observations. For Abell 2163 and Abell 2218, although the isothermal sphere with a finite 
core cannot produce giant arc images, a density model with a central cusp can produce a 
hnite Einstein radius, which is smaller than the observed radii. We find that a total mass 
density profile which declines as ~ produces the largest radius in models which are 
consistent with the X-ray surface brightness prohle. As the result, the extremely large ratio 
of the lensing mass to the X-ray mass is improved from 2.2 to 1.4 for Abell 2163, and from 3 
to 2.4 for Abell 2218. Eor RX J1347.5-1145, which is a cooling flow cluster, we cannot reduce 
the mass discrepancy. 

Subject headings: dark matter — galaxies:clusters;general — gravitational lensing — X-rays: 
galaxies 



1. Introduction 


Masses of clusters of galaxies has been estimated based on observations of X-ray gas and 
by observation of gravitational lensing. Although cluster mass estimated by these methods 
should yield the same values, they generally are are not, that is, the lensing mass is 2 ~ 3 
times larger than the X-ray mass(Grossman & Narayan 1989). The X-ray mass reflects the 
gas distribntion weakly bonnd by the gravitational force. The predicted mass distribution 
has a large core radius. Conversely, the lensing mass distribution is strongly concentrated to 
produce distorted images around the center. 

Many attempts have been made to reconcile the discrepancies between the X-ray and 
lensing mass measnrements have been made. Mass models with ellipticity can rednce the 
mass reqnired by the gravitational lensing (Miralda-Escnde 1993). However, the decrease is 
not so signihcant, especially, for Abell 1689 and Abell 2218 (Miralda-Escude & Babul 1995). 
Miralda-Escude & Babul (1995) presented a comparative analysis of lensing masses and X-ray 
masses, and found that a temperature two times higher than the observed value used be 
required to reconcile the discrepancy. Makino (1996) investigated the temperature distribution 
of the gas in lensing clusters, and showed that the gas temperatnre should signihcantly increase 
toward the center. The pressnre of the diffuse magnetic held could support a gravitational 
force strong enough to produce giant arcs (Loeb & Mao 1994). However, this reqnirement for 
the magnetic held leads to higher Faraday rotation measures in the lensing cluster than in 
nearby clusters (Makino 1997). Moreover, Loeb & Mao (1994) proposed that in addition to 
the gas pressnre, the bulk how or turbulent motions in the cluster structure would contribute 
to the support of the X-ray gas. Gas dynamical simulations, however, cannot prodnce strong 
dynamical motions (Evrard, Metzler, & Navarro 1996). 

Wu & Fang (1996) pointed out, by comparing the masses for samples of clusters derived 
from the X-ray to masses for samples from the gravitational lensing clusters, that the mass 
discrepancy is more signihcant aronnd the cluster center than at large radii. Indeed, the 
discrepancy between the X-ray mass and the weak lensing mass is smaller than that between 
X-ray mass and the strong lensing mass which is derived from giant arc images. Moreover, the 
discrepancy in cooling how clusters, which show strong peaks in the X-ray surface brightness 
at their centers, can be resolved by taking acconnt of the ehect of the peak in the snrface 
brightness on the mass estimate (Allen, Fabian, & Kneib 1996; Allen 1997). This analysis is 
not applicable to the non-cooling how clusters. 

The cluster density distribution has so far been believed to be smooth and hat at the 
cluster center. Indeed, both isothermal spheres with hnite core radii p oc l/(r^ -|- r^) and 
modihed King models, p oc l/(r^ -|- have been used to estimate the cluster mass. 

However, using N-body simulations of clusters of galaxies, Navarro, Frenk, & White (1997) 
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showed that density cusps form at the center, and that the density prohles have the same 
shape, independent of halo mass and of initial density fluctuation spectrum. According to 
them, the central density prohle is well £t by an analytic form, p oc l/[r(r + rs)^] (hereafter 
referred to as the NFW profile). The NFW profile is found in N-body/hydrodynamic 
simulations (Navarro, Frenk, & White 1995; Eke, Navarro, & Frenk 1997). Furthermore, 
Carlberg et al. (1997) found from measurements of galaxy velocity dispersions in clusters that 
the NFW profile can reproduce the measured velocity dispersion. At smaller scales, Hernquist 
(1990) derived an analytically approximated density profile, p oc l/[r(r + Vg)^] from the de 
Vaucouleurs law for the surface brightness of ellipticals. The Hernquist profile is in good 
agreement with numerical simulations (Dubinski & Carlberg 1991). 

However, the value of the power of the density cusp at the center is controversial because of 
the affect of gravitational softening and poor resolution at small scales. Fukushige & Makino 
(1997) showed from simulations with particle numbers one order of magnitude higher, that 
typical central density profiles are shallower than p ~ but steeper than p ~ r~^. Moore 
et al. (1997) also derived steeper inner density profiles, p ~ from cosmological N-body 
simulation with higher resolution. Furthermore, Evans & Collett (1997) found a steady-state, 
self-consistent, cusped solution of the collisional Boltzmann equation corresponding to 
p ~ Their solution is not thermodynamically stable but is dynamically stable. 

In this paper, we report on a study of the effect of density cusps in clusters of galaxies 
on the lensing mass and the X-ray mass. Assuming simple cluster mass models with cusps, 
i.e. spherically symmetric clusters with isothermal gas, we explicitly show that it is possible 
that the density cusp at the cluster center reduces the mass discrepancy between the X-ray 
mass and the lensing mass. In §2 we show our two step method: first we fix the parameters 
in the mass density models by comparing the estimated X-ray surface brightness profiles 
with the observed ones, and second, derive the estimated Einstein radii in the models and 
compare them with the observed Einstein radii from arcs. In §3, our method is applied to the 
observation of the lensing clusters, Abell 2163, Abell 2218 and RX J1347.5-1145, the X-ray 
masses and lensing masses are derived and the mass discrepancy are examined. In §4 we 
summarize our conclusions, and discuss our results. 

Throughout this paper, we assume Hq = 50kms“^ Mpc“^, Hq = 1 and Aq = 0. 


2. Method of analysis 

We consider a spherically symmetric cluster of galaxies with an X-ray emitting gas cloud, 
which is massive and distant enough to produce giant arcs around the center. We assume that 
the gas cloud of temperature T is isothermal and in hydrostatic equilibrium, and that the 
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total mass density profile p{r) in the clnster can be described by a cnsped density profile, 


p{r) 


do 

x^(l + 


( 1 ) 


where po is the normalization parameter and x = r/rp, Vg is the scaling parameter. For 
example, the mass density with (p, z/. A) = (0,2,1) is an isothermal sphere with a hnite core 
radins and (1,1,2) corresponds to the NFW prohle. N-body simnlations show that mass 
densities at the clnster scale decline as r~^ at large radii (Navarro, Frenk, Sz White 1995; 
Navarro, Frenk, & White 1996). We adopt the parameters that satisfy p + uX = 3, except for 
an isothermal sphere. 

The gas density distribntion, pg, satishes the equation of hydrostatic equilibrium: 

kT dlnpg GM{r) 

3 9 ’ \^) 

/Xm^p ^ 


where pm and rup denote the mean molecular weight (we adopt pm = 0.59 below) and the 
proton mass, respectively. Here we have neglected the contributions of gas and galaxies to 
the total mass M(r) within a radius r. The total mass M(r) is estimated by integrating 
equation (|^), and the gas density prohle is uniquely determined from equation (j^), if the gas 
temperature is given by observation. The gas density is expressed by two parameters, rg and 
h, dehned as 


27kT 


( 3 ) 


For the NFW prohle (p = 1, iz 
analytic form: 

Pgir) 


1, A = 2) one hnds that the gas density prohle has the 


Pgo exp 



In (1 + a;)\ 

X j ’ 


( 4 ) 


where pgo is the central gas density (Makino, Sasaki, & Suto 1997). 

To compare directly our gas density models with observation, we compute the X-ray 
surface brightness given by 


I{e) oc 





rdr, 


( 5 ) 


where Di is the angular diameter distance of the cluster from the observer. We have introduced 
a hnite cut-oh radius rcut to avoid divergence in equation (^. Since a mass density prohle 
which falls oh as r~^ at large radii cannot bind the isothermal gas, the gas density is hat at 
large radii. The hat gas density at inhnity leads to divergence in the surface brightness. 

A simpler mass density prohle has been extensively used to obtain the cluster mass from 
X-ray observations (e.g., Henry, Briel, & Nulsen 1993; Elbaz, Arnaud, & Bohringer 1995). 
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The observed X-ray emission is conventionally fitted with an isothermal /3-model (see, Sarazin 
1988) 


I{9) = Jo 


1 + 


' e 

Jr. 


1 


( 6 ) 


where Jq is the central surface brightness, 9c is the angular core radius of the X-ray surface 
brightness, and (3 is the slope parameter. In the /3-model, the gas density given by from 
equation (|^) with r^ut = oo is: 

- 3 / 3/2 


Pg(r) = pgo 


1+1 — 

Tr 


( 7 ) 


where Tc = 9cD\. Substituting equation into equation (H), one obtains the total mass. 
Furthermore, differentiating equation with respect to r, one finds the total mass density 
prohle 

3/3A:T 3 + {r/rcf 


p{r) = 


( 8 ) 


AnGpcarripr^ [1 -t- (r/rc)‘^f 
This density distribution is different from the cusped density distribution in equation (|^. The 
mass density in equation (^) is non-singular at the center. The cluster with a mass density 
given by (||) would not have a large enough mass to produce giant arcs around the center. 

On the other hand, we know that clusters do produce giant arcs. This imposes constraints 
on possible mass distributions, independent of X-ray mass estimates. In spherically 
gravitational lensing, if a lens is on the line of sight to a background source, the source will 
appear as arc images around a radius 9e at which the flux of a lensed image is apparently 
divergent in the limit of the geometrical optics. We refer to the radius 9^ as the Einstein 
radius. In observations of gravitational lensing we assume that the Einstein radius, 9^, is 
equal to the distance from the lens center to the giant arcs. 

Now for a cusped density profile we determine the parameters b and 9s in equation (|I]) and 
(§) on the basis of X-ray data on the surface brightness, and we estimate the cluster mass and 
the X-ray Einstein radius, 6 'e,x as follows. The X-ray estimated Einstein radius is independent 
of the observed radius 6*e, and the difference between the two reflects the mass discrepancy. 
In the case of a spherically symmetric lens, the Einstein radius is completely described by the 
surface mass density S(6*) dehned by 

p{r) 


fOO 

m = 2 

JD,e 


and the critical surface mass density 


Err — 


r2 - {Di9y 


Ds 


-rdr. 


( 9 ) 


4irG B|Bi. 


( 10 ) 
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where and Di^ are the angular diameter distances from the observer to the source and from 
the lens to the source, respectively. We dehne the projected mass within 6 by the parameters 
b and Vg determined from X-ray observation 


/ 

Jo 


■e 


m{e) = 27r / T.{e')e'de'. 


( 11 ) 


The X-ray estimated Einstein radius should satisfy 




( 12 ) 


where we use p{r) with the parameters determined from the X-ray observation. 

We dehne the projected lensing mass within the measured Einstein radius using the 
observed quantities, Scr, 0-^, and Di, as follows: 


ini.n,(SE) = 


(13) 


Note that the dehnition of the projected lensing mass is applicable to the region 6 *e because 
the giant arcs rehect the mass distribution within the measured Einstein radius only. Thus, 
the projected X-ray mass compared to the projected lensing mass is only constrained by the 
region within the measured Einstein radius 9-^. 

We can estimate an Einstein radius from the density models. We dehne the projected 
X-ray mass by the mass surface density E, determined from the X-ray observation: 



(14) 


where we again use p(r) with the above parameters determined by the X-ray observations. 
This is the X-ray mass to be compared with the lensing mass. 

For the NEW prohle, we can analytically derive m(0E) from the mass distribution. The 
detailed derivation of the Einstein radius for the NEW prohle was found in Bartelmann (1996) 
and Maoz et ah (1997). However, for most of our density models, the projected mass cannot 
be analytically estimated from the mass distribution. We numerically integrate the right hand 
of equation (0) to estimate the Einstein radius 6'e,x- 

Our method is summarized as follows: we compute the X-ray surface brightness for the 
various density mass models described by equations (2) and (|[). The estimated surface 
brightness is htted to the X-ray observation of lensing clusters to hx the parameters in the 
density models. Once the parameters are determined, we can compare the measured Einstein 
radius with the Einstein radius estimated from X-ray observation. We also estimate the 
projected mass and the Einstein radius for each density model. 


3. Results 


To study the effect of central mass concentration on the X-ray surface brightness and 
gravitational leasing, we adopt the following parameters in equation (|l]): v = 1 group; 
(/i,j/,A) = (1,1,2), (1.2,1,1.8), (1.4,1,1.6), (1.6,1,1.4), (1.8,1,1.2), and A = 1 group; (1,2,1), 
(1.2,1.8,1), (1.4,1.6,1), (1.6,1.4,1), (1.8,1.2,1). In addition to the cusped mass density, we 
compute the projected mass and the Einstein radius from the mass density derived from 
the /3-mo del in equation (^, and from the isothermal sphere with a finite core radius for 
comparison. The cut-off radius is fixed as rcut = lOrg in our calculation of equation (^). For 
the isothermal model (p, i/. A) = (0,2,1), we exceptionally fix r^t as lO'^rg. The gas densities 
derived from the cusped density profiles are flat at radii > lOrg. Since the cusped density 
profiles decline as the gas density cannot be strongly bound. Unless the gas distribution 
at large radii declines faster than corresponding io (3 = 0.4 in the /3-model, the gas 
density at the radius r ~ D\d mainly contributes to the integral in equation (^. As will be 
shown below, the cut-off radius of lOrg is larger than the radius in which X-ray emissions from 
the clusters are detected. Thus, the choice, r^ut = lOrg, is reasonable for our aim. The general 
behavior of the surface brightness estimated from the cusped density profile for the v = 1 
group is discussed in Suto, Sasaki, & Makino (1998). 

We apply our method of estimations to three clusters, Abell 2163, Abell 2218 and RX 
J1347.5-1145. They are X-ray clusters in which giant arcs have been found. Their cluster 
redshifts, gas temperature and observed Einstein radii are summarized in Table 1. We also 
list the parameters of the /3-model, /3 and dc-, derived from our fitting results of the observed 
surface brightness in Table 1. 


3.1. The X-ray surface brightness 

We illustrate the behavior of the surface brightness derived from cusped density profiles. 
In Figure 1 we show the X-ray surface brightness estimated from the cusped profiles for 
b = 0.7. The surface brightness for density models with z/ = 1 is plotted in Figure 1(a). 
The solid line indicates the surface brightness for (/r, z/. A) = (1,1, 2). The dotted line, short 
dotted line, long dotted line and dashed dotted line show the surface brightness profiles for 
(/r, u, A) = (1.2,1,1.8), (1.4,1,1.6), (1.6,1,1.4), and (1.8,1,1.2), respectively 

The surface brightness profiles with the parameters /i = 1, 1.2 and 1.4 are similar in 
shape. They show a flat core around the center and a steep power law shape at large radii. 
Makino et al. (1998) showed that the gas density for /r = 1 is in good agreement with the 
/3-model in equation (|^) which gives a good fit to the surface brightness of real clusters of 
galaxies. Thus, the surface brightness profiles for /i = 1, 1.2 and 1.4 are expected to agree 
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with the observation of clusters of galaxies. However, the surface brightness for /i = 1.8 shows 
a single power law shape and no hnite core. It would be difficult to £t the surface brightness 
for fi = 1.8 to observations. For fi = 1.6 the surface brightness is marginal. The core around 
the center deviates slightly from a flat prohle. While its surface brightness could agree with 
cooling flow clusters, it is not a good fit to clusters with a flat and large core radius like the 
Coma cluster. 

The decrease of the core radius for larger fi is due to the strong concentration of the mass 
density at the center for large /x. Especially for /i = 1.8 the core radius is not found in Figure 
1(a). The behavior of the surface brightness for fi = 1.8 is substantially the same as /x = 2. 
The gas density in the density model with /i = 2 is divergent at the center. In order to see the 
behavior of gas density for /i = 2, we expand the integrand in equation (|[) around the center, 
X <C 1: 


dlnpg dvrCpohm^p ^ 


dlnr 


kT 


^(1 — \x’^)dx. 


( 16 ) 


If /i = 2, the right hand in equation (ll^) is hnite at the limit r —0. The gradient of the gas 
density is hnite at the center. The density model with p = 2 cannot have a hat core in surface 
brightness under the assumption of the isothermal gas. For p = 1.6, the X-ray core radius is 
extremely small. 

In Figure 1(b) we plot the surface brightness for the parameters (/i, v, A) = (1.2,1.8,1), 

(1.4.1.6.1) , (1.6,1.4,1), and (1.8,1.2,1). The surface brightness in the A = 1 model is steeper 
than in the v = 1 model at large radii for the same value of h. It is evident from equation 
(p!5|) that a cusped density prohle with a smaller u and a larger A has a larger core radius. 
For example, the surface brightness for (/i, i/, A)= (1,1,2) in Figure 1(a) is hatter than that for 

(1.2.1) in Figure 1(b). 


3.2. Abell 2163 

Abell 2163 is a distant and elliptical cluster of galaxies. It is classihed as an Abell richness 
class 2 cluster.The redshift is z\ = 0.201. Abell 2163 has a very high temperature of 12.4keV 
(Holzapfel et ah 1997) and X-ray luminosity in the 2-10 keV band of Lx = 6.0 x 10^® ergs s“^ 
(Elbaz et ah 1995). Two arcs have been observed around the central galaxy. The measured 
Einstein radius is 15."6, and its redshift is 0.728 (Miralda-Escude & Babul 1995). 

We ht the surface brightness of the cusped density profiles to the observations by the 
minimization, and determined the values of b and 6s- In Figure 2, we plot the observed X-ray 
surface brightness prohles of Abell 2163 derived from the ROSAT HRl images and the best-ht 
prohles for various {p, z/. A) and the /3-model. 

The center of the surface brightness prohle is located on the X-ray maximum of the cluster 
emission. The errors are 1-cr errors and account for Poisson statistics. A hat exposure map is 
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assumed. That is, vignetting effects are not taken into account in the data. The original data 
analysis can be found in Elbaz et ah (1995). 

All the models are good hts at large radii, but the behavior of the prohle for each of our 
models is very different in the flat core region. The contribution of the density cusp to the 
surface brightness prohle is signihcant in the hat core region. Thus, we conhne our attention 
to the behavior of the best-ht prohles within below. 

The best-ht prohles for the = (0,2,1), (1,1,2) and (1.2,1,1.8) models are very 

similar each other. They approximately reproduce the core structure of the X-ray emission. 
These three prohles deviate slightly upwards from the /5-model within 6c toward the center. 
The surface prohle for (/i, z/. A) = (1.4,1,1.6) model is steeper than the models with /i < 1.2, 
but is within the error bars. The X-ray prohles of more centrally concentrated models, 

(/U, z/. A) = (1.6,1,1.4) and (1.8,1,1.2) are too steep in the central region to agree with the 
innermost data point within the error bars. The best-ht prohles of these models become 
steeper as fi increases. This tendency is also seen in Figure 2(c) for the A = 1 group. The 
behaviors of the prohles for the A = 1 group within dc are slightly hatter than ones for the 
z/ = 1 group for the same value of /U. The deviation from the best-ht /5-model for the A = 1 
group is relatively small. 

From the reduced values in Table 2(a) and Figure 3(a), we hnd that the /5-model is the 
best htting model for the X-ray emission in Abell 2163, and that the cusped density prohles 
for fi = 1, 1.2 and 1.4 marginally agree with the data. It is apparent from Figure 3(a) that the 
reduced values for the cusped models change over the range 1.3 ~ 1.8, except for fi = 1.6 
and fi = 1.8. Although we present the reduced only for the z/ = 1 model, the behavior of 
the reduced x^ fhe A = 1 group is similar. The relatively large values of the reduced x^ ^ire 
caused by substructures in Abell 2163 which appear in the X-ray image (Elbaz et ah 1995). 
If we could get rid of the contribution of the substructures, the models for /i < 1.4 would be 
in agreement with the data. 

We estimate the X-ray predicted Einstein radii 6 'e,x and the projected masses mx within 
the measured Einstein radius Oe for Abell 2163. The projected lensing mass of Abell 2163 is 
estimated to be 4.13 x IO^^Mq. The results are summarized in Table 2. For the z/ = 1 model 
the mass ratio and Einstein radius ratio are plotted in Figure 3 (b) and (c). For the z/=l 
group, the ratios of the projected lensing mass to the projected X-ray mass monotonically 
approach one as fi increases, but none is less than one. Although a larger fi decreases the 
mass ratio m\ens/mx, the deviation of the estimated X-ray surface brightness aronnd the 
center from the data is signihcantly large. This tendency is also found in the ratio ^E,X/^E. 
The ratio increases as fi increases for the z/ = 1 gronp. Interestingly, a hnite Einstein radius 
can form in the cusped density prohles althongh the giant arcs cannot be produced in the 
density models with a hnite core. This implies that the mass concentration around the center 
is reqnired to produce giant arcs. From the X-ray and the lensing analysis, onr best mass 
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density model is (/i, u, A) = (1.4,l,1.6) for Abell 2163. This density profile increases toward the 
center as p oc The estimated Einstein radius for the our best model amounts to 47% 

of the observed value. The projected mass for our best model is 57% larger than the mass 
estimated from the best-£t /3-model. 

For the A=1 group, we obtain the same results as for the u = 1 group from Table 2(a). 
We, however, note that the change in b for z/=l groups is less sensitive to /i than in the A=1 
groups. The mass discrepancy is not improved for the A=1 group as significantly as for the 
z/=l group. 


3.3. Abell 2218 

Abell 2218 is an Abell richness class 4 cluster at 2 : 1 = 0 .175. The radial velocity dispersion 
of cluster galaxies is very high, 1370kms“^ (Le Borgne, Pello, & Sanahuja 1992). About 30 
arcs have been identified in this cluster. The configuration of these arcs strongly suggests 
that this cluster has two mass clumps (Kneib et al. 1995). While there is a secondary bright 
clump around galaxy ^244 in a optical observation (Pello et al. 1992), the X-ray image shows 
only one maximum close to the position of the cD galaxy. We consider one giant arc ^(^359 
with a radius of 20."8. The redshift of the source galaxy is 0.702 (Miralda-Escude & Babul 
1995). This arc is mainly caused by the central clump around the cD galaxy. Thus, we can 
neglect the second mass clump. The X-ray luminosity. Lx, is 0.7 x 10"^® ergs s“^ (0.1-2.4 
keV), and the gas temperature is measured to be 7.2keV (Mushotzky & Loewenstein 1997). 
Birkinshaw & Hughes (1994) reported 6c = 60" and (3 = 0.65. On the other hand, Markevitch 
(1997) reported that the X-ray emission within a core radius of 60" is resolved into several 
components. Markevitch speculated that this may be caused by lensed X-ray emission, a 
merger shock or a gas trail of an infalling subgroup, or both lensing and a merger. Under 
these assumptions, he reanalyzed the X-ray profile and obtained 9c = 26" and (3 = 0.49. We 
do not consider this possibility. We simply fit the surface brightness data with our models 
and we ignore the substructures in the cluster. 

In Figure 4, we show the observed X-ray surface brightness prohles of Abell 2218 derived 
from the ROSAT HRI images and the best-fit profiles for various values of (/i, v, A). The 
analysis of the data for Abell 2218 is the same as for Abell 2163. The detailed description of 
the data analysis is found in Squires et al. (1996) 

We also plot the best-£t /3-model with f3 = 0.56 and 6c = 53". The data points are 
smoothly distributed in comparison with the data of Abell 2163. The behavior of the best-fit 
prohles of the v = 1 and A = 1 groups for Abell 2218 is basically similar to those for Abell 
2163. All the profiles are in excellent agreement with the data in the power law region at 
large radii. All of the surface brightnesses estimated from the cusped models are steeper than 
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that from the best-ht /9-model at the flat core region. The best-£t profile becomes steeper as 
jj, increases. All of the profiles of the u = 1 group deviate from the error bars. The surface 
brightness of the A = 1 group pass inside the error bars only for /i < 1.2. 

We list results of the lensing analysis for Abell 2218 in Table 2(b). For the u = 1 model, 
the results are plotted in Figure 3. The projected lensing mass of Abell 2218 is 6.3 x lO^^M©. 
The mass ratio and the ratio of the Einstein radius change for the parameters b and Vg in 
similar fashion as for Abell 2163. The change of b and 6s for Abell 2218 shows the same 
behavior in Table 2(a) for both u = 1 and A = 1 groups. The reduced values are smaller 
than the values for Abell 2163 because of the smoothness of the X-ray emission. As seen 
in the mass estimate for Abell 2163, both the isothermal sphere and the best fit-/9-model 
mass models cannot produce giant arcs for Abell 2218. The cusped density models, however, 
produce finite Einstein radii, however these are less than the estimated radius. For cusped 
density models, the ratio, miens/'^rx and 6 *e,x/6*e approach one as /i increases. For Abell 2218, 
our best model is (/i, z/, A) = (1.4,l,1.6). The estimated Einstein radius for our best model 
for Abell 2218 corresponds to 16% of the observed value. The projected mass mx for our 
best model is 27% larger than the /9-model mass. We cannot reduce the mass discrepancy of 
Abell 2218 as significantly as we did for Abell 2163. This is mainly due to the different gas 
temperature of the two clusters. 


3.4. RX J1347.5-1145 

A distant {z\ = 0.451) cluster RX J1347.5-1145 has two dominant galaxies. It is the most 
X-ray luminous cluster: Lx = 7.3 x 10^^ ergs s“^ (0.1-2.4 keV) and T=9.3 keV(Schindler et 
ah 1997). The X-ray profile of this cluster has a strong peak. Schindler et ah (1997) noted 
that this strong peak is probably evidence of cooling flows. However, they could not find 
significant evidence for cool gas in the spectra. Two bright arcs are located at about 35" from 
the cluster center and are positioned on opposite sides of the central galaxy (Schindler et ah 
1995). The redshift of these arcs is 0.81 (Sahu et ah 1998). 

In Figure 5, we present the observed surface profile of RX J1347.5-1145 derived from the 
ROSAT WRl images (Schindler et al. 1997) and the best-ht profiles estimated from various 
models. As a result of fitting by the /9-model, we find an extremely small core radius Vc 
of 56kpc in comparison with that of nearby clusters, ~ 250kpc (Jones & Forman 1984). 
Indeed, the core radii of Abell 2163 and Abell 2218 yield about 200-300kpc in our htting. The 
observed profile is steeper than the the best-ht /9-model inside the core radius. This steep 
prohle with a small core radius is presumably evidence of cooling hows, as Schindler et al. 
(1997) suggested. As shown in Figure 5, all of the best-ht profiles computed from the cusped 
density profiles is steeper than the best-ht /9-model. The surface brightness for /x < 1.6 of 
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both z/ = 1 and A=1 group are in excellent agreement with the data, because of the steepness 
of the data points at the center of this cluster. 

The projected lensing mass of RX J1347.5-1145 is 6.26 x The uncertainty of 

the mass is caused by the unknown redshift of the lensed source galaxy. Table 2(c) indicates 
that the reduced values in RX J1347.5-1145 are approximately equal to one for both the 
/i = 1 and u = 1 groups, except for the models with (i = 1.8. Thus, models with jj, < 1.6 
are acceptable. However, the mass discrepancy is not decreased for RX J1347.5-1145. The 
estimated Einstein radius and the projected X-ray mass are almost unchanged, even if /i 
increases. The projected mass mx for (/i, z/, A) = (1.6,l,1.4) model is only 7% larger than the 
/3-model. This may be due to the small core radius of this cluster. The X-ray core radius 9^ of 
this cluster is smaller than the Einstein radius 0 e, although the X-ray core radius is large than 
the Einstein radius for the previous two clusters. As shown in Figure 5, the surface brightness 
prohles are almost unchanged at 6* > 9^ even if the model is changed. Therefore, we conclude 
that the projected mass discrepancy of RX J1347.5-1145 was not decreased by our method 
based on available ROSAT data. The ROSAT HKl data may be dominated by X-rays from 
the cooling flow, not from the cluster as a whole, however. To improve our analysis for RX 
J1347.5-114, we need the X-ray data from the cluster as a whole, and we should take cooling 
flows into account. 


4. Conclusions and Discussion 

We studied the effect of the density cusp on the mass cluster estimates. We compared the 
Einstein radii estimated from the density prohle with a cusp with the observed radius. We 
showed that the marginally steep cusp, /i = 1.4, reduces the discrepancy of the mass estimates 
for Abell 2163 and Abell 2218. The ratio of the lensing mass to the X-ray mass is improved 
from 2.2 to 1.4 for Abell 2163, and from 3 to 2.4 for Abell 2218. 

We found that the clusters should have a density cusp at the center to reduce the 
discrepancy between the X-ray mass and the lensing mass. The cusped density model is 
consistent with the density prohle found in the N-body simulations (Navarro et al. 1997). 
Navarro et al. (1997) concluded that the density prohle behaves as ~ r~^ toward the center. 
However, the more strongly cusped density prohle at the center is expected to reconcile the 
mass discrepancy. In fact, our best model falls oh p oc r~^'^ around the center. The density 
prohles of clusters of galaxies might be more strongly peaked at the center as Fukushige & 
Makino (1997) and Moore et al. (1997) have shown with N-body simulations. 

Although the introduction of a cusp in the density prohle is not sufficient to resolve the 
mass discrepancy, the density cusp mainly contributes to the reduction of the discrepancy. In 
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addition to density cusps, we should take into account the gas temperature decline, the second 
mass clump and the gas dynamical motion. We expect that ellipticity in the mass distribution 
plays an important role in resolving the discrepancy. Indeed, Bartelmann (1996) has shown 
that mass estimation under the assumption of a spherically symmetric mass distribution 
are likely to overestimate the mass by ~ 30%-50% . If the density cusp model reduces the 
ratio of the lensing mass to the X-ray mass miens/'^rx to 1.3 — 1.6, the mass discrepancy 
would be resolved by the adoption of a more detailed cluster model, especially elliptical mass 
distributions for non-cooling flow clusters (Asano & Makino 1998). For RX J1347.5-1145, the 
difference between the lensing mass and the X-ray mass was not improved. This may be due 
to cooling flows (Schindler et al. 1997). For this cluster the cooling flows should be taken into 
account (Allen 1997). 

Furthermore, independent of X-ray observations it has been recognized that the cusped 
models with 1 < /i < 1.4 are more plausible in gravitational lensing analysis. Bartelmann 
(1996) pointed out that the NFW prohle predicts large values of the radial magnihcation for 
the cluster of galaxies MS 2137. This problem implies that the corresponding sources must be 
surprisingly thin in the radial direction. Evans & Wilkinson (1997) have shown that lenses 
which have more singular density cusps than the NFW prohle assuage Bartelmann’s problem. 
The result of Evans & Wilkinson, which was derived from the analysis of gravitational lensing, 
agrees well with our result. 

We adopted a density model which declines as ~ r~^ at large radii. However, the density 
prohle at large radii is uncertain. In N-body simulations it is difficult to determine whether 
mass clumps at large radii are assigned to the cluster. In fact, the prohle of the Hernquist 
model for elliptical galaxies decreases as ~ r~^ at large radii. Dubinski & Carlberg (1991) 
have shown by numerical simulations that the prohle of clustering objects falls oh as ~ r~^ 
at large radii. Furthermore, Markevitch & Vikhlinin (1997) have applied the mass density 
model p oc a;“^(l -|-a;)^““ to Abell 2256, and ht the data with the surface brightness estimated 
from mass models for a < 5. However, mass densities which decline faster than r~^ are less 
ehective for the formation of giant arcs than our model. 

Miralda-Escude & Babul (1995) have studied how cusped density models resolve the mass 
discrepancy. They have shown that a cusped mass density produces a peak in the isothermal 
gas density around the center, and that the cusped density model is inconsistent with 
ROFATPSPC data of Abell 2218, Abell 2163 and Abell 1689. The ROFATPSPC data show 
extremely small errors. On the other hand, we have used the ROSAT HRl data, which have 
relatively large errors. The two detectors are very different in the spatial resolution and their 
backgrounds. In addition, the ht is sensitive to the existence of substructures in the clusters, 
and the results depend on the data. If the surface brightness estimated from our models is ht 
to the RO^ATPSPC data, the mass discrepancy would not be signihcantly reduced. 

Lastly, we comment on the scaling parameter in the NFW prohle. By numerical 
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simulations of dark halos, Navarro, Frenk, & White (1996) predicted that the scaling 
parameter Vg is approximately several hundred kpc at the cluster mass scale. On the other 
hand, Vg derived from our analysis are about IMpc for the Abell 2163 and 2218 clusters. 
Namely, the size of Vc predicted by Navarro, Frenk, & White is smaller than the observed 
values. This problem has been already pointed out by Makino et al. (1997). Both Bartelmann 
(1997) and Evans & Wilkinson (1997) adopted the scaling parameter which satishes the 
estimation by the numerical simulations of Navarro, Frenk, & White. Their results for MS 
2137 will be altered by taking into account the constraint by X-ray observations. 


We are grateful to Makoto Hattori, Yasushi Suto and Kenji Tomita for their helpful 
advices, and Doris Neumann for providing the X-ray surface brightness data of Abell 2163 
and Abell 2218. We thank Sabine Schindler for kindly providing the ROSAT WRl data of RX 
J1347.5-1145 in ASCII format. Lastly, we appreciate Maxim Markevitch for kindly providing 
the ROSAT B.RI images for Abell 2218. 
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Figure captions 

Fig. 1.— The X-ray surface brightness prohles estimated from the mass density models. The 
parameter b is equal to 0.7. (a)the surface brightness prohles for the u = 1 group. The solid 
line indicates the surface brightness for the models with the parameters (/i, z/, A) = (1,1,2). 
The dotted line, the short dashed line, the long dashed line, and the dotted-dashed line show 
the surface brightness for the models (/r, u, A) = (1.2,1,1.8), (1.4,1,1.6), (1.6,1,1.4), (1.8,1,1.2), 
respectively, (b) the surface brightness prohles for A = 1 group. The solid line indicates the 
surface brightness for the models with the parameters (/i, z/. A) = (1, 2,1). The dotted line, the 
short dashed line, the long dashed line, and the dotted-dashed line show the surface brightness 
for the models (/i, z/, A) = (1.2,1.8,l), (1.4,1.6,1), (1.6,1.4,1), (1.8,1.2,1), respectively. 

Fig. 2.— The X-ray surface brightness prohle of the X-ray emission from Abell 2163. The boxes 
are the data of the observed surface brightness. The ROSAT HRl data were kindly provided 
to us by Doris Neumann. The solid lines are the best-£t /5-model with parameters, (3 = 0.54 
and dc = 65". (a) We present the surface brightness for the isothermal sphere with a finite core 
radius model and the z/ = 1 group for /x < 1.2. The dotted line indicates the surface brightness 
for the isothermal sphere (/i, z/. A) = (0, 2,1). The short dashed line and the long dashed line 
show the surface brightness for the cusped models with parameters (/i, z/. A) = (1,1,2) and 
(1.2,1,1.8), respectively, (b) The X-ray surface brightness for the z/ = 1 groups with /i > 1.2. 
The dotted line, short dashed line and the long dashed line show the surface brightness for the 
cusped models with parameters (/i, z/. A) = (1.4,1,1.6), (1.6,1,1.4), and (1.8,1,1.2), respectively, 
(c) The X-ray surface brightness for the A = 1 groups with /i < 1.4. The dotted line, short 
dashed line and the long dashed line show the snrface brightness for the cusped models with 
parameters (/r, z/. A) = (1,2,1), (1.2,1.8,1), and (1.4,1,6.1), respectively. 

Fig. 3.— The rednced valnes, the ratios of the X-ray estimated Einstein radius to the 
observed Einstein radius and the ratios of the lensing mass to the X-ray mass for Abell 2163 
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and Abell 2218. We present the results only for the v = 1 group. The behavior for the A = 1 
group is similar to the z/ = 1 group. The solid lines and dotted lines indicate the quantities 
of Abell 2163 and Abell 2218, respectively, (a) the reduced values for various values of the 
parameter p. (b) the ratios of the X-ray estimated Einstein radius to the observed Einstein 
radius, (c)the ratios of the lensing mass to the X-ray mass 

Fig. 4.— The X-ray surface brightness for Abell 2218. The ROSAT HRI data were kindly 
provided to us by Doris Neumann. The parameters for the best-£t /5-model are (3 = 0.56 
and 9c = 53". We show the parameters for the other models in Table 2(b). (a) The surface 
brightness for the isothermal sphere with a finite core radius model and the v = 1 group for 
/r < 1.4. The dotted line indicates the surface brightness for the isothermal sphere. The short 
dashed line, the long dashed line and the short and dashed line indicate the surface brightness 
profiles for the cusped models with parameters (/r, z/. A) = (1,1,2), (1.2,1,1.8), and (1.4,1,1.6) 
respectively, (b) The X-ray surface brightness for the A = 1 groups with fx < 1.6. The dotted 
line, short dashed line, the long dashed line and the short and long dashed line show the surface 
brightness for the cusped models with parameters (/i, z/. A) = (1,2,1), (1.2,1.8,1), (1.4,1,6.1) 
and (1.6,1.4,1), respectively. 

Fig. 5.— The X-ray surface brightness for RX J1347.5-1145. The ROSATBBJ data were kindly 
provided to us by Sabine Schindler. The parameters for the best-fit /5-model are f3 = 0.56 and 
dc = 8.2". In Table 2(c), we show the parameters for the other models, (a) The surface 
brightness for the isothermal sphere with a finite core radius model and the v = 1 group for 
fi < 1.4. The dotted line indicates the surface brightness for the isothermal sphere. The short 
dashed line, the long dashed line and the short and dashed line indicate the surface brightness 
profiles for the cusped models with parameters (/i, z/. A) = (1,1,2), (1.2,1,1.8), and (1.4,1,1.6) 
respectively, (b) The X-ray surface brightness in RX J1347.5-1145 for the A = 1 groups with 
fi < 1.6. The dotted line, short dashed line, the long dashed line and the short and long dashed 



line show the surface brightness for the cusped models with parameters = (1,2,1), 

(1.2,1.8,1), (1.4,1,6.1) and (1.6,1.4,1), respectively. 



Cluster 

2:1 


T 

9-e 

Oc 

/9 




keV 

arcsec 

arcsec 


Abell 2163 

0.201 

0.728 

12.4 

15.6 

65 

0.54 

Abell 2218 

0.175 

0.702 

7.2 

20.8 

53 

0.56 

RX J1347.5-1145 

0.451 

0.81 

9.3 

35 

8.2 

0.56 


Table 1. Cluster sample. The parameters of (3 and 9c are derived from our analysis. 


(/i, z/, A) 

h 

9c/9s 

9e,x/9e 

^lens/ 

reduced y' 

/5-model mass 

- 

- 

- 

2.24 

1.32 

(0,2,1) 

0.15 

1.9 

- 

2.01 

1.55 

(1,1,2) 

0.63 

0.25 

0.10 

1.75 

1.51 

(1.2,1,1.8) 

0.56 

0.15 

0.27 

1.62 

1.55 

(1.4,1,1.6) 

0.54 

0.064 

0.47 

1.43 

1.76 

(1.6,1,1.4) 

0.49 

0.015 

0.72 

1.23 

2.63 

(1.8,1,1.2) 

0.39 

0.00043 

0.96 

1.04 

5.02 

(1,2,1) 

0.31 

0.37 

0.033 

1.95 

1.37 

(1.2,1.8,1) 

0.32 

0.25 

0.18 

1.75 

1.44 

(1.4,1.6,1) 

0.35 

0.12 

0.42 

1.50 

1.69 

(1.6,1.4,1) 

0.36 

0.032 

0.70 

1.24 

2.58 

(1.8,1.2,1) 

0.34 

0.00067 

0.89 

1.10 

5.12 


Table 2(a). Results for Abell 2163 
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(/i, z/, A) 

b 

Oc/ds 

6*E,x/6*E 

"^lens/ 

reduced 

/9-model mass 

- 

- 

- 

3.02 

0.59 

(0,2,1) 

0.15 

1.8 

- 

2.78 

0.66 

(1,1,2) 

0.66 

0.24 

0.013 

2.69 

0.65 

(1.2,1,1.8) 

0.58 

0.15 

0.070 

2.58 

0.67 

(1.4,1,1.6) 

0.57 

0.060 

0.16 

2.38 

0.76 

(1.6,1,1.4) 

0.52 

0.014 

0.29 

2.18 

1.26 

(1.8,1,1.2) 

0.39 

0.00047 

0.42 

2.00 

2.76 

(1.2,1) 

0.33 

0.37 

0.0025 

2.93 

0.61 

(1.2,1.8,1) 

0.33 

0.25 

0.045 

2.74 

0.63 

(1.4,1.6,1) 

0.37 

0.11 

0.14 

2.47 

0.73 

(1.6,1.4,1) 

0.37 

0.030 

0.28 

2.20 

1.24 

(1.8,1.2,1) 

0.35 

0.00063 

0.39 

2.11 

2.82 


Table 2(b). Results for Abell 2218 
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(/i, z/, A) 

h 

Oc/ds 


miens /mx 

reduced 

/9-mo del mass 

- 

- 

0.27 

2.89 

1.03 

(0,2,1) 

0.14 

2.3 

0.27 

2.82 

1.00 

(1.1.2) 

0.60 

0.29 

0.25 

2.96 

1.00 

(1.2,1,1.8) 

0.53 

0.19 

0.24 

2.89 

1.00 

(1.4,1,1.6) 

0.45 

0.098 

0.25 

2.81 

1.01 

(1.6,1,1.4) 

0.41 

0.028 

0.26 

2.70 

1.09 

(1.8,1,1.2) 

0.58 

0.000089 

0.25 

3.00 

1.66 

(1.2.1) 

0.32 

0.39 

0.23 

3.12 

1.02 

(1.2,1.8,1) 

0.31 

0.27 

0.24 

2.98 

1.01 

(1.4,1.6,1) 

0.31 

0.16 

0.24 

2.87 

1.01 

(1.6,1.4,1) 

0.30 

0.056 

0.25 

2.77 

1.08 

(1.8,1.2,1) 

0.47 

0.00031 

0.26 

2.93 

1.65 


Table 2(c). Results for RX J1347.5-1145 
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